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Abstract 

We present an alternative way to determine the unknown pa- 
rameter associated to a gaussian approximation in a generic two- 
dimensional model. Instead of the standard variational approach, 
we propose a procedure based on a quantitative prediction of con- 
formal invariance, valid for systems in the scaling regime, away from 
criticality. We illustrate our idea by considering, as an example, the 
sine-Gordon model. Our method gives a good approximation for the 
soliton mass as function of [5. 



Pacs: ll.lO.Lm 

ll.lO.Kk 



" Institute de Fisica La Plata, Departamento de Fisica, Facultad de Ciencias Exactas, 
Universidad Nacional de La Plata. CC 67, 1900 La Plata, Argentina, 
emails: iucci@fisica.unlp.edu.ar, naon@lisica.unlp.edu.ar 



1 



The so called "self consistent harmonic approximation" (SCHA) is a non- 
perturbative technique that has been extensively employed in Statistical Me- 
chanics [0 and Condensed Matter physics P |§] applications. Roughly 
speaking it amounts to replacing an exact action Stmc by a trial action S'trial 
that makes the problem tractable. Usually Striai is just a quadratic action 
that depends on certain unknown parameter Q that must be determined 
through some criterion such as the minimization of the free energy of the 
system. This approximation is intimately related the the "gaussian effective 
potential" (GEP) in Quantum Field Theories (QFT's), a variational 

approximation to the effective potential which uses a gaussian wave func- 
tional depending on some mass parameter as the trial ground state. It also 
relies on a minimization principle often called "principle of minimal sensi- 
tivity" [§] to determine the additional parameter. The main purpose of this 
Letter is to point out that in two-dimensional problems (1+1 QFT's) there 
is an alternative way to obtain the quantity fl. Since, as we shall see, this 
method is based on Conformal Field Theory (CFT) [0, we call it "confor- 
mal gaussian approximation" (CGA). In view of the interesting physics that 
is currently under investigation in low- dimensional systems (organic conduc- 
tors [11|, charge transfer salts |T^, quantum wires |]13[, edge states in a 



two-dimensional (2D) electron system in the fractional quantum Hall regime 
H|, Carbon Nanotubes (CNT) [fH], Luttinger-like systems in 2D high tem- 



perature superconductors [0]) the dimensional restriction in the validity of 
our approach does not diminish its practical interest. 

We shall begin by depicting the main features of the standard SCHA. 
One starts from a partition function 



•'true 



where T>fi is a generic integration measure. An elementary manipulation 
leads to 
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for any trial action Striai- Now, by means of the property 
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for / real, and taking natural logarithm in equation (Q), we obtain Feynman's 
inequality 
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In general, S'trial depends on some parameters, which are fixed by minimizing 
the right hand side of the last equation. 

From now on, in order to illustrate the procedure, we shall consider scalar 
2D models, with action 



St 



{2tt 



(5) 



where (p{p) is a scalar field and F{p) is usually of the form F{p) ~ p^. For 
simplicity, in this formula we have written the kinetic term in Fourier space 
but we kept the interaction term in coordinate space. 
As the trial action one proposes a quadratic one. 
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where Vt is the trial parameter. To proceed with the minimization, we first 
write the partition function for the trial action in the form 



-Atrial = exp 



V 
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where V is the volume (infinite) of the whole space J (Px, and we have defined 
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with the formal properties 
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To go further we must specify the potential V{ip). Let us specialize the 
discussion to the sine-Gordon model (SGM), taking 



a 



v{^) = ^ [1 - cosm] 



(12) 



and F{p) = p^. 

As it stands the theory is divergent. To take into account the divergences 



we can implement Coleman's normal order prescription [|T8| with respect to 
any given mass-dimensional constant p from the beginning [|l^, by replacing 
in the lagrangian V{ip) by its normal ordered form 



a 
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where A/'p[...] means the normal ordering form with respect to p. 

It is now straightforward to compute (Struc — •S'triai), by following, for 
instance, the steps explained in ref. PO]. The result is 



(5'truo 'Atrial) trial ~ ^ 
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Now inserting equations (|^ and ([T^) in equation (^, and extremizing the 
r.h.s. with respect to fi, we finally obtain 
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^2 _^^-PV2{hm-h{p)) (15) 

This gap equation allows to extract a finite answer for fi, depending on the 
mass parameter p (the difference Ii{Vt) — Ii{p) is finite). Note that the value 
of p is completely arbitrary, if one chooses it to be equal to the trial mass fl, 
the solution to the equation is 



n"^ = a. (16) 
The same result would be obtained if instead of p = f2 one had taken p = y/a. 



Let us now present an alternative route to determine Q. To this end we 
will exploit a quantitative prediction of conformal invariance for 2D systems 
in the scaling regime, away from the critical point. Indeed, starting from the 
so called 'c-theorem' Cardy ||2^ showed that the value of the conformal 



anomaly c, which characterizes the model at the critical point, and the second 
moment of the energy-density correlator in the scaling regime of the non- 
critical theory are related by 



rf^x|xP(.(xH0)> = 3^^jp-^. (17) 

where e is the energy- density operator, is its scaling dimension and t oc 
(T — Tc) is the coupling constant of the interaction term that takes the 
system away from criticality. The validity of this formula has been explicitly 
verified for several models [H il H [H- For the SGM, the energy density 



operator is given by the cosine term, its conformal dimension is = /3^/47r, 
t is the coupling constant a/P"^ and the associated free bosonic CFT has 
c= 1. 

Now we claim that Q can be determined in a completely different, not 
variational way, by enforcing the validity of the above conformal identity for 
the trial action. In other words, we will demand that the following equation 
holds: 



— / £x (cos/5<^(a;) cos pip {0)),^-^^ = - — (18) 
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which is to be viewed as an equation for the mass parameter f2. Of course, 
if one is interested in comparing the answer given by this formula with the 
SCHA result, when evaluating the left hand side of (p^Sf ) one must adopt a 
regularizing prescription equivalent to the normal ordering implemented in 
the SCHA calculation. A careful computation leads to the following gap 
equation: 



p 32 

where we have defined the variable u = jS'^/Air {0 < u < 2) and p is the 
normal ordering parameter, as before. We see that, as in the standard SCHA 
equation (plSf) , one has different answers for different choices of p, but in this 
case, the results obtained for the values y/a and Q are different. In any case 
one gets a non trivial dependence of fl on jS"^ in contrast with the SCHA. This 
is interesting if one recalls the physical meaning of mass gaps in the context 
of the SGM. Indeed, as it is well-known, Dashen, Hasslacher and Neveu 
(DHN) have computed by semiclassical techniques the mass spectrum 
for the SGM. It consists of a soliton (associated to the fermion of the Thirring 
model) with mass 



Msoi = v^, (20) 

TT U 

and a sequence of doublet bound states with masses 
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with N = 1,2,... < (2 — u)/u. (From this last condition it is easy to see 
that in order to have N bound states one must have u < 2/{N +1). As a 
consequence there is no bound state for u > 1). Therefore the masses in the 
SGM spectrum also depend on u as the CGA mass of equation (|T9|). Thus, in 
this respect our proposal seems to be able to improve the standard gaussian 
prediction for the SGM, at least qualitatively. In order to perform a more 
specific and quantitative discussion let us compare equations (0) and (pOD 
as functions of u. We set p = ^Ja^ which corresponds to the prescription 
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employed by DHN when deriving (^) and (PT]). The result is shown in Fig. 
1 where one can observe a general qualitative analogy between both curves. 
In particular, for 0.7<m<1(m = 1 corresponds to the free fermion point 
of the Thirring model and to the Luther-Emery point in the backscattering 
model |2^) our CGA prediction is in full agreement with the values of the 
soliton mass as computed by DHN. We want to stress that for m = 1 we get 
Vl/^/a = a/3/32 0.30 whereas the exact value given by ( pOD is l/vr ^ 0.31 
(standard SCHA yields, of course, Vt/\/a = 1). 



To conclude, we have reconsidered the well-known SCHA method in which 
a comparatively complex action is replaced by a simpler quadratic system 
depending on a mass parameter Q which is usually determined through a 
variational calculation. Taking into account the (l+l)-dimensional case, we 
have proposed an alternative way for evaluating Q. Our technique (CGA) 



is based on a consequence of Zamolodchikov's |2T| c-theorem first derived 
by Cardy [^. We have illustrated the proposal by considering the sine- 
Gordon (SGM) model. We showed that for this model CGA gives a quite 
good prediction for the behavior of the soliton mass as function of (see 
equations ( [I9| ) and ( pOD and Fig. 1. It would be interesting to test our 
approach in other models such as the continuum version of the tricritical 
Ising model, which is described by the second model of the unitary minimal 
series 



28[] [ p9[] with central charge c = 7/10. 
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Figure caption 



Figure 1: Masses in units of ^Ja as functions of u. The dashed hne is 
Mgoil \/oL^ whereas the filled line represents Vt/\/a as given by CGA. 
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